PROGRAM MAPS
USE POLYMORPHIC_COMPLEXTAYLOR
TYPE(DAMAP) M1,M2,M2oM1
integer mf
mf=20
open(unit=mf,file='results.txt')
CALL INIT(NO1=2,ND1=1,NP1=0,NDPT1 =0)     !   <------------------ init for maps in ND1 degrees of freedom

CALL ALLOC(M1,M2,M2oM1)





M1%V(1)= 1.d-1 + (2.D0.MONO.'10') - (3.D0.MONO.'01') + (1.D0.MONO.'20') 
M1%V(2)= 2.d-1 + (1.D0.MONO.'10') + (1.D0.MONO.'01') 

M2%V(1)= 2.d-1 + (2.D0.MONO.'10') 
M2%V(2)=         (1.D0.MONO.'01') 

M2oM1=M2.o.M1


write(20,*) " Map M1 "
CALL PRINT(M1,mf)
write(20,*) " Map M2 "
CALL PRINT(M2,mf)

write(20,*) " Map M2oM1  "
CALL PRINT(M2oM1,mf)



CALL KILL(M1,M2,M2oM1)

close(20)
END PROGRAM MAPS

